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ABSTRACT 


This  report  contains  an  analysis  of  the  antenna  calibration 
data  taken  with  a  high  frequency  (HF)  antenna  array  and  a  pilot 
source  towed  behind  a  small  aircraft.  The  analysis  used  an  array 
calibration  method  developed  by  E.K.L.  Hung.  The  objective  was  to 
determine  whether  this  method  was  adequate  for  calibrating  the  HF 
antenna  arrays  in  the  Canadian  Forces  Supplementary  Radio  (CFSRS) 
modernization  project. 

The  study  showed  that  the  calibration  method  significantly 
improved  the  performance  of  the  array  in  detecting  low- elevation 
(11°  or  less)  targets  and  estimating  their  directions.  The  degree 
of  improvement  decreased  as  the  target  elevation  increased.  At  an 
elevation  of  49°,  it  was  marginal  at  best. 

This  study  also  uncovered  the  presence  of  systematic  errors  in 
the  calibration  data.  These  errors  increased  rapidly  with  the  pilot 
source  elevation. 

It  is  not  possible  to  decide  whether  or  not  the  calibration 
method  is  adequate  for  the  CFSRS  antenna  arrays.  In  this  study,  the 
calibration  method  could  be  inadequate  for  high-elevation  targets, 
or  the  method  was  adequate  but  the  systematic  errors  negated  the 
improvements  gained  in  the  calibration. 


RESUME 


Ce  rapport  presents  une  analyse  des  donnees  de  calibration 
recueillies  avec  un  reseau  d'antennes  decametriques  et  une  source 
temoin  remorquee  derriere  un  petit  avion.  L' analyse  utilise  une 
methods  de  calibration  d'antennes  developpee  par  E.  K.  L.  Hung. 
L'objectif  de  cette  analyse  est  de  determiner  si  cette  methods  est 
adequate  pour  la  calibration  d'antennes  decametriques  utilisees 
dans  le  pro jet  de  modernisation  du  reseau  radio  supplementaire  des 
Forces  canadiennes  (RRSFC) . 

L' etude  demontre  que  la  methods  de  calibration  ameliore 
considerablement  la  performance  de  I'antenne  lors  de  la  detection 
de  cibles  a  petit  angle  de  site  (11°  ou  moins)  et  lors  de 
1' estimation  de  leurs  directions.  Le  facteur  d  amelioration 
decroit  lorsque  1' angle  de  site  de  la  cible  augments.  A  un  angle 
de  site  de  49°,  1 ' amelioration  est  tout  au  plus  marginals. 

Cette  etude  a  aussi  permis  de  decouvrir  la  presence  d'erreurs 
systematiques  dans  les  donnees  de  calibration.  Ces  erreurs 
augmentent  rapidement  avec  1' angle  de  site  de  la  source  temoin. 

II  n'est  pas  possible  de  determiner  si  la  methods  suggeree  est 
adequate  pour  le  reseau  d'antennes  decametriques  utilise  dans  le 
pro jet  de  modernisation  du  RRSFC.  Dans  cette  etude,  la  methode  de 
calibration  est  soit  inadequate  pour  de  grands  angles  de  site  ou 
elle  est  adequate  mais  les  erreurs  systematiques  annulent  les 
ameliorations  obtenues  dans  la  calibration. 
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EXECUTIVE  SUMMARY 


This  report  contains  an  analysis  of  the  calibration  data  taken 
with  a  high  frequency  (HF)  antenna  array  in  1994  at  Leitrim,  near 
the  Ottawa  airport.  It  uses  an  array  calibration  method  developed 
by  Eric  K.L.  Hung.  The  objective  is  to  determine  whether  this 
method  is  adequate  for  calibrating  the  HF  antenna  arrays  in  the 
Canadian  Forces  Supplementary  Radio  Systems  (CFSRS)  modernization 
project.  These  arrays  are  used  for  detecting  long-range  HF  signals 
and  estimating  their  directions. 

Calibration  optimizes  the  performance  of  antenna  arrays  in 
target  detection  and  direction  finding. 

The  method  developed  by  Hung  is  unique.  It  is  highly  tolerant 
to  the  presence  of  interfering  signals,  including  unknown  specular 
multipath  signals.  Other  calibration  methods  do  not  have  this 
property.  They  cannot  be  used  to  calibrate  the  CFSRS  arrays, 
because  interfering  signals  are  always  present  in  HF  antenna 
environments . 

The  following  observations  are  based  on  the  analysis  of  data 
files  generated  with  transmitter  ranges  not  less  than  7  km  and 
transmitter  elevations  not  higher  than  11°. 

1.  With  a  few  exceptions,  the  calibration  method  increases  the 

signal-to-interference-plus-noise  ratio  (SINR)  in  the  array 
snapshots  by  a  factor  of  2  to  4  (3  to  6  dB) .  Because  of  this 

general  increase,  the  method  helps  us  detect  signals  that  are 
undetectable  before  calibration.  It  also  reduces  the 
integration  time  needed  to  raise  the  SINRs  of  the  desired 
signals  to  the  minimum  detectable  level . 

2.  The  calibration  method  increases  the  accuracy  of  target 
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direction  estimates.  More  specifically,  it  reduces  the  standard 
deviation  of  the  errors  in  bearing  estimates  by  a  factor  of  3.0  or 
more . 

3  .  The  method  increases  the  resolution  performance  of  the  antenna 
array.  On  the  average,  it  reduces  the  number  of  array 
snapshots  with  unresolvable  equal-strength  signals  by  a  factor 
of  2  or  more . 

Other  data  files,  generated  with  shorter  transmitter  ranges 
and  higher  transmitter  elevations,  were  available  for  the  analysis. 
They  required  corrections  for  the  perturbations  in  element  response 
patterns  as  well  as  for  the  difference  in  transmitter  elevations  at 
the  array  elements.  These  corrections  could  not  be  carried  out. 
Consequently,  the  authors  hesitate  to  draw  conclusions  on  the 
results  they  produced,  although  these  results  generally  agreed  with 
the  above  observations. 

It  is  not  yet  possible  to  decide  whether  the  calibration 
method  developed  by  Hung  is  adequate  for  the  HF  antenna  arrays  in 
the  CFSRS  modernization  project.  The  decision  requires  the  analysis 
of  calibration  data  generated  with  transmitter  elevations  higher 
than  11°  and  transmitter  ranges  not  shorter  than  7  km.  These  data 
are  currently  not  available. 

The  difficulties  in  using  the  higher  elevation  files  prompted 
the  authors  to  study  an  alternative  calibration  method  that 
tolerates  the  dependence  of  calibration  coefficients  on  directions. 
Preliminary  results  indicated  that  this  alternative  method  was 
superior  to  the  method  in  this  report.  Because  of  these  results, 
the  authors  have  decided  to  stop  further  work  on  the  method  in  this 
report  and  direct  the  remaining  efforts  to  the  evaluation  of  the 
alternative  method. 
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1 . 0  INTRODUCTION 


An  accuirat©  calibration  of  an  antonna  array  is  ©ssontial,  if 
th©  array  is  to  b©  us©<i  in  <iir©ct ion -finding'  with  high- r© solution 
(HR)  methods  such  as  th©  MUSIC  method  [1]  .  Otherwise,  th© 
performance  of  these  HR  methods  in  estimating  and  resolving  signal 
directions  will  be  degraded. 

Many  methods  to  calibrate  a  narrowband  antenna  array  have  been 
published  in  the  open  literature.  They  include  those  by  Paulraj  and 
Kailath  [2],  Friedlander  and  Weiss  [3,4],  Brown,  McClellan  and 
Holder  [5,6],  and  Pierre  and  Kaveh  [7].  However,  these  methods 
cannot  be  used  in  array  environments  where  unknown  interfering 
signals,  including  diffuse  and  specular  multipath  signals,  are 
present.  In  particular,  they  cannot  be  used  to  calibrate  high- 
frequency  (HF)  antenna  arrays,  because  multipath  interfering 
signals  are  always  present. 

A  unique  calibration  method  [8]  has  recently  been  developed 
for  HF  antenna  arrays  at  the  Defence  Research  Establishment  Ottawa 
(DREO) .  The  method  equalizes  the  gains  and  phases  of  the  elements 
in  antenna  arrays.  Unlike  the  methods  in  [2]  to  [7]  ,  it  is  highly 
tolerant  to  the  presence  of  interfering  signals,  including  unknown 
diffuse  and  specular  multipath  signals.  In  preliminary  studies, 
this  method  produced  highly  accurate  estimates  of  the  element  gains 
and  phases. 

This  report  contains  an  analysis  of  the  calibration  data 
measured  in  1994  with  an  HF  antenna  array  at  Leitrim,  near  the 
Ottawa  airport.  It  uses  the  calibration  method  in  [8]  .  The  main 
objective  is  to  determine  whether  the  method  is  adequate  for 
calibrating  the  HF  antenna  arrays  in  the  Canadian  Forces 
Supplementary  Radio  Systems  (CFSRS)  modernization  project. 
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The  remainder  of  this  report  is  organized  as  follows:  Chapter 
2  contains  a  description  of  the  antenna  array  and  the  coordinate 
systems  used  in  the  analysis.  Chapter  3  describes  the  calibration 
measurements  and  the  data  collected.  Chapter  4  discusses  the 
preprocessing  of  these  calibration  data.  It  includes  checking  the 
possible  presence  of  faulty  elements,  Doppler  filtering,  and 
correcting  the  deviation  of  signal  wavefronts  from  spherical 
wavefronts.  The  calibration  method  is  presented  in  Chapter  5,  and 
the  estimates  of  element  gains  and  phases  obtained  with  this  method 
are  in  Chapter  6.  Chapter  7  provides  an  evaluation  of  the 
calibration  method;  it  examines  the  effect  of  calibration  on  the 
accuracy  of  direction  estimates,  the  signal-to-interference-plus- 
noise  ratio  (SINR)  in  the  array  snapshots,  and  the  resolution  of 
signal  directions.  The  conclusions  are  given  in  Chapter  8. 
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2.0  ANTENNA  ARRAY 


The  antenna  array  was  located  at  Leitrim,  which  is  near  the 
Ottawa  International  Airport.  It  was  planar,  had  twelve  elements, 
and  was  part  of  a  larger  array  with  52  elements.  Each  one  of  these 
elements  was  a  six-meter  high  vertical  monopole,  implemented  by  a 
metal  pipe  standing  on  eight  ground  radials  and  fed  approximately 
one-third  of  the  way  up.  The  positions  of  these  elements  are  given 
in  Table  1  and  are  plotted  in  Figure  1. 

There  are  two  coordinate  systems  in  Table  1:  a  UVW- rectangular 
system  supplied  with  the  data  and  a  polar  system  to  be  used  in  the 
analysis  of  calibration  data.  These  systems  are  defined  in  Figure 
2 .  In  the  rectangular  system,  the  U-direction  is  the  north 
direction,  the  V-direction  is  west,  and  the  W-direction  is  up.  In 
the  polar  system,  the  location  of  a  point  is  identified  by  a 
triplet  {r,d,\p).  Here,  r  is  the  distance  of  the  point  from  the 
origin,  6  is  the  bearing  angle,  and  ^  is  the  elevation  angle.  The 
bearing  6  is  measured  in  the  UV-plane  and  clockwise  from  the  U- 
axis . 


The  rectangular  coordinates  (u,v,w)  and  polar  coordinates 
(r,d,\l/)  of  a  point  in  space  are  related  to  each  other  by: 

u  =  r  cosi/'  COS0  , 

V  =  -  r  cos^  sin0  , 
w  =  r  sintp  , 

and 

r  =  (u^+v^+w^)  , 

6  =  tan‘^(-v/u)  , 

xj/  =  tan'^  (w/ (u^+v^)  . 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 


The  bearing  d  is  related  to  the  azimuth  angle 


in  an  XYZ- 
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coordinate  system,  defined  with  {X, Y,Z}= {East, North, Up} ,  by 


=  90°  -  e  .  (2.7) 

The  coordinates  in  the  UVW  and  XYZ  systems  are  related  by 

X  =  -V  ,  (2.8) 

y  =  u  ,  (2.9) 

z  =  w  .  (2.10) 


ELEMENT  POSITIONS 

N 

Rectangle  Coordinate 

Polar  Coordinate 

u 

V 

w 

r 

e 

1 

3.263 

-24.786 

o 

o 

25 

82.5 

0 . 0 

2 

-3.263 

-24.786 

0.0 

25 

91 .5 

0.0 

3 

-15.219 

-19.834 

o 

o 

25 

127.5 

o 

O 

4 

-19.834 

-15.219 

o 

o 

25 

142.5 

O 

O 

5 

-24.786 

3.263 

o 

o 

25 

187.5 

0.0 

6 

-19 . 834 

15.219 

o 

o 

25 

217.5 

O 

O 

7 

23 . 097 

9.567 

o 

o 

25 

337.5 

0.0 

8 

24 . 786 

3.263 

o 

o 

25 

352.5 

0.0 

9 

85.001 

22.776 

o 

o 

88 

345.0 

O 

O 

10 

-22.776 

85.001 

o 

o 

88 

255.0 

0.0 

11 

120 . 741 

32.352 

0.0 

125 

345.0 

0.0 

12 

-32.352 

120.741 

o 

o 

125 

255.0 

O 

O 

Table  1.  Element  positions  in  the  HF  antenna  array. 
Distances  u,  v,  w,  and  r  are  measured  in  meters.  Angles 
d  and  rp  are  measured  in  degrees . 
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V(west) 


ISO 


Figure  1 .  Geometry- 
array  . 
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U(north) 


3 . 0  CALIBRATION  MEASUREMENTS 


In  the  calibration  measurements,  a  vertically  oriented  dipole 
transmitter  was  towed  behind  an  airplane.  This  plane  was  flown 
along  specified  paths  while  snapshots  of  the  outputs  of  the  array 
elements  were  taken.  Two  global  positioning  systems  were  used  to 
determine  the  ranges  and  directions  of  the  transmitter  in  the 
snapshots.  The  first,  identified  as  GPSl,  was  at  the  antenna  site. 
Its  output  time,  denoted  by  tl,  was  recorded  with  the  array 
snapshots.  The  second,  identified  as  GPS2,  was  a  differential 
system  consisting  of  two  GPS  receivers,  one  on  the  plane  and  the 
other  on  a  known  ground  location.  Its  output  time  t2  and  position 
rGps2(t2)  were  recorded  in  a  log  file  on  the  aircraft  at 
approximately  one  second  intervals . 

Six  transmitter  frequencies  were  used:  5.1,  7.5,  9.3,  11.5, 
15.1,  and  18.0  MHz.  The  outputs  at  each  array  element  were 
downshifted  in  frequency  to  near-baseband  and  sampled  with  an  A/D 
converter  at  10.4  kHz.  Approximately  fifty  data  files  were 
generated.  Table  2  contains  the  nominal  flight  parameters  and  sizes 
for  the  files  generated  at  5.1  MHz.  Each  data  block  has  8000  array 
snapshots.  Figures  3  to  10  illustrate  the  paths  of  the  airplane  in 
these  files.  They  include  five  clockwise  (CW)  circles  round  the 
array,  one  counter-clockwise  (CCW)  circle  round  the  array,  and  two 
straight  lines  over  the  array. 

The  flight  path  in  File  2  has  a  side  loop  near  bearing  210° 
and  a  hook  near  bearing  333°.  Their  presence  are  not  in  the  flight 
plan.  In  later  discussions,  we  shall  refer  to  File  2  as  the 
modified  file  with  the  side  loop  and  the  hook  removed. 

The  transmitter  was  towed  behind  the  airplane.  Its  location  at 
the  recorded  GPS2  time  t2  was  approximately  the  same  as  that  of 


6 


GPS2  at  time  t2-A,  where  A  was  equal  to  the  GPS2-to-transmitter 
distance  divided  by  the  speed  of  the  airplane.  The  value  of  A  was 
approximately  2  seconds;  its  exact  value  depended  on  the  wind 
speed,  wind  direction,  plane  speed,  and  plane  direction  at  time  t2  . 

In  a  separate  study  of  the  calibration  data  by  Dr.  R.W  Jenkins 
[9] ,  it  was  noted  that  the  recorded  GPSl  time  tl  was  not  equal  to 
the  GPS2  time  t2,  recorded  at  the  same  time.  These  times  were 
related  by 


where  e  was  a  non- zero  difference  independent  of  frequency  and 
t3^ansmitter  elevation.  The  origin  of  this  difference  was  unknown. 
This  study  also  showed  that  the  transmitter  position  at  time  tl, 
denoted  by  rxx(tl),  could  be  derived  from  the  GPS2  position  at  t2 
as 


rTx(tl)  —  ^GPS2 

=  r3ps2(tl-A')  ,  (3.2) 

where  A'=e+A.  The  assignment 

A'  =  8.7  seconds  (3.3) 

produced  the  best  overall  agreement  between  the  estimated  and 
measured  transmitter  directions. 
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File 

Date 

Elevation 

(deg) 

Plannee 

Altitude 

(ft) 

i  Track 

Radius 

(km) 

Direction 

No .  of 

Blocks 

1 

18/05/94 

4.5 

3000 

10 

CW 

2255 

2 

20/05/94 

11 

5000 

7 

CW 

2308 

3 

20/05/94 

11 

5000 

7 

CCW 

1438 

4 

24/05/94 

18 

6000 

5 

CW 

1106 

5 

24/05/94 

31 

6150 

3 

CW 

723 

6 

24/05/94 

49 

6000 

2 

CW 

367 

7 

24/05/94 

6000 

1350 

8 

24/05/94 

6000 

1365 

Table  2 .  Nominal  flight  parameters  for  the  data  files 
measured  at  5.1  MHz.  The  date  format  is  day/month/year . 
Each  data  block  contains  8000  array  snapshots 
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Figure  3.  Flight  path  in  Fil4  1.  The  airplane  circles 
the  array  with  a  radius  of  10  km  and  proceeds  clockwise 
(CW)  through  405°,  from  north  45°  east  to  north  90  east. 


U(km) 


the  side  loop  at  bearing  210°  and  hook  at  bearing  at 
333°.  The  airplane  circles  the  array  with  a  radius  7  km, 
proceeds  CW  through  358°,  from  north  24°  west  to  north 


U(k: 


Figure  5 . 
(CCW)  405° 


Flight  path  in  File  3 
,  north  19°  west  to  nc 


U(kr 


Figure  6 .  Flight  path  in  File 
south  49°  west  to  north  57°  west 


0 


1 


i;  counter-CW 
it . 


start 


Figure  10.  Flight  path  in  File  8.  It  is  a  straight  1 
over  the  array  with  a  heading  of  north  1.0®  east. 


4 . 0  DATA  PREPROCESSING 


The  data  files  were  preprocessed  before  they  were  used  to 
calculate  the  element  gains  and  phases.  The  following  is  a  list  of 
the  tasks  in  this  preprocessing: 

1.  Check  if  there  are  faulty  elements  in  the  generation  of  the 
data  files; 

2 .  Reduce  the  number  of  array  snapshots  in  each  file  by  Doppler 
filtering;  and 

3 .  Correct  the  deviations  of  signal  wavefronts  from  planar 
wavefronts  in  the  filtered  snapshots. 

Task  2  is  required,  because  the  files  are  extremely  large  and  have 
2.9  to  19  million  array  snapshots  per  file. 

Three  corrections  are  ignored  in  the  preprocessing.  The  first 
is  the  correction  for  the  difference  in  the  orientations  of  the 
array  elements.  It  has  been  assumed  that  the  metallic  pipes  used  in 
the  array  elements  were  vertical  and,  therefore,  parallel.  The 
relative  orientations  of  some  pipes  were  checked,  using  the 
parallax  method,  in  the  summer  of  1995.  They  were  not  exactly 
parallel,  although  the  largest  angle  between  any  two  pipes  was  at 
most  a  few  degrees.  This  correction  is  ignored,  because  the 
orientations  of  the  elements  were  not  measured  in  the  calibration 
experiments.  The  second  correction  is  for  the  direction-dependent 
differences  in  the  gain  and  phase  response  patterns  of  the  array 
elements.  These  patterns  were  not  identical  in  their  direction 
dependence  for  the  calibration  measurements,  although  all  were 
approximately  dipole.  This  correction  is  ignored,  because  it 
depends  on  many  unknown  factors  such  as  the  ground  conductivity. 
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vegetation,  local  structures,  etc.  The  third  correction  is  for  the 
difference  in  the  transmitter  elevations  at  the  array  elements.  The 
transmitter  range  was  10  km  or  less  in  the  calibration  experiment. 
Consequently,  its  elevation  at  each  of  the  array  elements  was 
different  in  any  particular  array  snapshot.  This  correction  cannot 
be  carried  out,  because  it  involves  the  unknown  element  patterns 
and  orientations. 

Ignoring  the  first  two  corrections  produces  perturbations  in 
the  gain  and  phase  response  patterns  of  the  array  elements. 
Usually,  these  perturbations  are  larger  if  the  transmitter 
elevation  is  higher.  Ignoring  the  last  correction  produces 
amplitude  and  phase  errors  at  the  outputs  of  the  array  elements. 
These  errors  are  larger  if  the  transmitter  range  is  shorter. 

File  6  in  Table  2  was  measured  with  nominal  transmitter  range 
and  elevation  of  2  km  and  49°,  respectively.  Because  the  range  was 
short  and  the  elevation  was  high,  we  do  not  expect  to  get  reliable 
results  from  this  file.  A  large  percentage  of  the  array  snapshots 
in  Files  7  and  8  was  measured  with  ranges  shorter  than  2  km  and 
elevations  higher  than  49° .  Consequently,  we  omit  these  files  in 
later  studies . 

4 . 1  Faulty  Element  Check 

The  calibration  measurements  were  such  that,  after  Doppler 
filtering  with  a  few  hundred  array  snapshots  or  more,  the 
transmitter  signal  at  the  array  elements  was  very  strong  compared 
with  all  other  signals.  Therefore,  the  power  spectrum  calculated 
with  the  outputs  of  each  array  element  should  have  a  strong 
amplitude  at  the  Doppler  shifted  frequency  of  the  transmitter 
signal .  It  follows  from  this  property  that  the  array  did  not  have 
faulty  elements  if  (a)  the  spectrum  of  each  element  had  a  very 
strong  signal  present  and  (b)  the  frequency  of  this  signal  was  the 
same  at  all  elements. 
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Some  results  obtained  in  a  check  for  faulty  elements  are 
presented  in  Figure  11.  The  input  data  used  in  Doppler  filtering 
are  the  first  1000  array  snapshots  in  Block  #1  of  the  first  5.1  MHz 
data  file.  There  are  twelve  spectra,  one  spectrum  from  the  outputs 
of  each  array  element.  The  strongest  amplitude  in  each  spectrum  is 
at  a  baseband  frequency  of  1.95  kHz  and  is  at  least  25  dB  stronger 
than  the  other  signals.  Therefore,  the  array  did  not  have  faulty 
elements  during  the  generation  of  this  5.1  MHz  file. 

Only  1000  array  snapshots  are  used  to  calculate  each  Doppler 
spectrum,  although  each  data  block  has  8000  snapshots.  The  choice 
on  the  number  of  snapshots  was  influenced  by  many  factors,  e.g., 
the  FFT  frequency  bin  width  must  not  be  smaller  than  the  bandwidth 
of  the  transmitter  signal,  the  transmitter  was  swinging  in  the  air, 
the  transmitter  antenna  orientation  was  changing  all  the  time,  etc. 

In  the  calculation  of  each  spectrum,  we  augment  the  1000  real 
output  amplitudes  of  an  element  with  1048  zeros,  convert  the 
augmented  sequence  to  a  2048  point  complex  frequency  spectrum  with 
the  Fast  Fourier  Transform  (FFT)  ,  and  plot  the  powers  of  the 
positive  frequency  points  only.  The  augmentation  with  zeros  is 
necessary.  Otherwise,  the  picket  fence  effect  and  noise  presence 
could  combine  to  prevent  the  strongest  amplitude  from  appearing  at 
the  same  frequency  bin  in  all  twelve  spectra. 


The  FFT  frequency  bins  are  identified  with  one  to  2048. 
Because  the  sampling  frequency  is  10.4  kHz,  the  frequency  bandwidth 
of  the  bins  is  10.4  Hz  (=  10.4  kHz  divided  by  1000,  the  number  of 
snapshots  used  in  Doppler  filtering)  and  the  separation  of  the  bins 
is  5.0781  Hz  (=  10.4  kHz  divided  by  2048,  the  number  of  FFT  bins) . 
Bins  1025  to  2048  are  negative  frequency  bins,  because  the  array 
snapshots  have  in-phase  (real)  components  only. 

The  above  1.95  kHz  peak  frequency  is  equal  to  Ifg-fol,  where 
is  the  frequency  of  the  Doppler  shifted  transmitter  signal  at 
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the  array  elements  and  is  the  frequency  of  the  local  oscillator 
in  down  conversion.  It  corresponds  to  the  FFT  frequency  bin  number 
385,  where  the  strongest  amplitude  is  present. 

There  are  three  easily  identifiable  interfering  signals  in 
each  spectrum.  Their  frequencies  are  0,  0.06,  and  3.08  kHz  and 
their  strengths  are  more  than  25  below  that  of  the  strongest  signal 
at  1.95  kHz.  The  0  Hz  interfering  signal  probably  comes  from  an 
offset  in  the  A/D  converters.  The  origins  of  the  others  are  not 
known . 

The  above  three  interfering  signals  are  also  present  in  all 
other  5.1  MHz  files.  An  example  is  given  in  Figure  12,  which  is 
calculated  with  the  first  1000  snapshots  in  Block  #1  of  the  second 
5.1  MHz  file. 
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Spectrum  (dB)  Spectrum  (dB)  Spectrum  (dB)  Spectrum  (dB) 


Figure  11.  The  spectra  for  each  of  the  twelve  array- 
elements,  using  the  first  1000  array  snapshots  in  Block 
1  of  File  1. 


Spectrum  (dB)  Spectaim  (dB)  Spectrum  (dB)  Spectrum  (dB) 


Figure  12 .  The  spectra  for  each  of  the  twelve  array 
elements,  using  the  first  1000  array  snapshots  in  Block 
1  of  File  2. 


4.2  Data  Reduction 


We  use  Doppler  filtering  to  reduce  the  size  of  each  data  file 
by  a  factor  1000.  The  procedure  described  below  takes  note  that  the 
files  are  stored  in  blocks  of  8000  array  snapshots.  There  are  two 
steps . 

Step  1. 

Separate  each  data  block  into  eight  sets  with  1000  consecutive 
snapshots  per  set. 

Step  2. 

For  each  set, 

(a)  Augment  the  1000  output  amplitudes  of  each  element  with 
1048  zeros; 

(b)  Calculate  the  2048  point  complex  FFT  frequency  spectra  of 
these  augmented  sequences; 

(c)  Study  the  spectrum  calculated  with  the  outputs  of  the 
first  element.  Identify  the  signal  bin  as  the  bin  with 
the  largest  power.  Ignore  negative -frequency  bins  1025  to 
2048;  and 

(d)  Construct  a  new  array  snapshot  as 

X  k  (Xi,X2,X3,  •  •  •  ,Xi2)'^  •  (4.1) 

In  (4.1),  Xi  to  Xi2  are  the  complex  FFT  amplitudes  in  the  largest- 
signal  bin  of  the  spectra  calculated  from  the  outputs  of  elements 
1  to  12,  respectively.  The  superscript  'T'  denotes  the  transpose. 
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The  above  procedure  preserves  the  relative  amplitudes  of  the 
transmitter  signal  at  the  array  elements.  It  also  maximizes  the 
signal-to-noise  ratio  (SNR)  in  the  array  snapshots.  Here,  the  ratio 
refers  to  the  noise  power  in  one  FFT  bin,  which  has  a  bandwidth  of 
10.4  Hz. 

In  later  discussions,  the  transmitter  location  for  x  will  be 
denoted  by  and  will  be  defined  as  the  average  of  the 
locations  in  the  1000  array  snapshots  used  to  construct  x. 

4.3  Correction  for  Non-Planar  Wavefront 

The  signal  wavefront  in  x  is  not  planar.  We  correct  the 
deviation  from  a  planar  wavefront  by  replacing  it  with  a  new 
snapshot  x  with  components  {x^, X2,  •  ■  • , X12}  given  by 

x„  =  (s„/rt)  {exp  [j27r(s„-rt) /X]  a„(0t,i/'t)}  ^ 

n=l,2, • ■ • ,12.  (4.2) 


In  this  expression, 

(si,  S2,  S3,  •  •  •  ,  S12}  =  element  distances  from  transmitter, 
X  =  signal  wavelength  (=  58.82  m  at  5.1  MHz)  , 

and 


a„(0t/(^t)  =  exp  [  j27r  (u„cosV'tCOs6t  "  v^cosi/'tsin^t 

+  WnSini/'t) /X]  .  (4.3) 

The  first  term,  (s^/r^)  ,  in  (4.2)  corrects  the  r'^  dependence  of 
signal  power  on  range;  the  second  term,  {•••},  equalizes  the  signal 
phases  at  the  elements;  and  a„(0t/^t)  is  the  nth  component  of  the 
array  steering  vector  for  direction  and  infinite  range. 
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5.0  CALIBRATION  METHOD 


A  brief  description  of  the  array  calibration  method  in  [8]  is 
presented  here.  It  contains  the  assumptions  made,  the  signal  model 
used,  the  theory  for  estimating  the  element  gains  and  phases,  and 
the  procedure. 

In  the  presentation,  the  set  of  element  gains  are  denoted  by 
the  vector  g  or,  equivalently,  the  set  {gn*  n=l,  2 ,  •  •  •  ,  12 }  . 
Similarly,  the  set  of  element  phases  are  denoted  by  <f)  or 
n=l,2,  '  •  •  ,12} .  For  convenience,  the  following  uniqueness  conditions 
are  imposed  on  the  gains  and  phases : 


gi  =  1  / 


(5.1) 


01  =  0  . 


(5.2) 


5.1  Assumptions 

The  following  assumptions  have  been  made  in  the  development  of 
the  calibration  method  in  [8]  : 

1.  The  relative  gain  responses  of  the  array  elements  are 
independent  of  elevation  and  azimuth, 

2 .  The  relative  phase  responses  of  the  array  elements  are 
independent  of  direction, 

3.  The  mutual  coupling  among  the  elements  can  be  ignored, 

4.  The  signal-to-interference-plus-noise  ratios  (SINRs)  at 
the  array  elements  are  reasonably  high. 
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5 .  The  calibration  signal  is  effectively  uncorrelated  with 
the  interfering  signals  in  the  data,  and 

6.  The  average  signal-plus-noise  powers  incident  upon  the 
array  elements  are  the  same. 

The  relative  gain  and  phase  responses  are  also  known  as  the  element 
gains  and  phases,  respectively. 

We  do  not  examine  the  validity  of  the  assumptions  in  this 
report.  We  only  study  what  one  can  gain  when  the  method  in  [8]  is 
used  to  calibrate  the  antenna  array  at  Leitrim. 

5.2  Signal  Model 

In  this  signal  model,  a  preprocessed  array  snapshot  has  the 
decomposition 

X  =  a  b  (g,  0,  ’At)  +  interference  +  noise.  (5.3) 

The  first  term  represents  the  direct  signal  from  the  transmitter. 
Here,  a  is  the  complex  signal  amplitude  at  a  reference  element  in 
the  array  and  b(g,0,0tf^t)  the  array  manifold  vector  for  the 

transmitter  direction  • 

The  array  manifold  vector  for  any  direction  (0,^)  is  given  by 

h(g,<l),  d  ,\p)  =  G(g)  T(0)  3l{6  ,  (5.4) 

where  G(g)  and  T (0)  are  NxN  diagonal  matrices  associated  with  the 
element  gains  and  phases,  respectively,  and,  under  the  assumptions, 
are  defined  as 


and 


G(g)  =  diag{gi,g2,  •  •  •  ,gi2} , 


(5.5) 
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T(0)  =  diag{exp  [j0i]  ,  exp  [j02]  /  •  ■  ■ '  }  • 


(5.6) 


5 . 3  Theory 


The  theoretical  developineiit  initially  discusses  the  estiination 
of  element  gains.  Next,  it  presents  a  transformation  that  equalizes 
the  gains  and  simultaneously  'fixes'  the  transmitter  direction  in 
the  data.  Finally,  it  explains  how  one  can  estimate  the  element 
phase  responses  with  the  transformed  data. 

5.3.1  Estimation  of  Element  Gains 

The  expected  values  of  the  incident  signal-plus-noise  powers 
are  equal  at  the  array  elements.  Consequently,  the  expected  values 
of  the  output  powers  at  these  elements  are  proportional  to  the 
squares  of  the  element  gains.  If  these  output  powers  are 
{P1/P2/ •  •  •  .P12}/  the  ratio  p„/g^  is  independent  of  the  element 
identification  n. 

Given  a  set  of  estimates  {Pn*  n=l ,  2 ,  •  •  •  ,  12 } ,  the  gain 
estimates  can  be  calculated  as 


,  n=l,2,---,12,  (5.7) 

J  Pi 

since  condition  (5.1)  requires  gi=l.  One  can  show  that  an  accurate 
estimation  of  the  average  output  powers  would  result  in  an  accurate 
estimation  of  the  element  gains. 

5.3.2  Transformation  for  Snapshots 


The  array  manifold  vector  for  the  transmitter  direction  is 
given  by  (5.4),  with  0  =  0t  and  Because  G(g)  and  T(</>)  are  both 
diagonal,  a  simple  rearrangement  leads  to 
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b{g,<^,  t(0). 


(5.8) 


where 


A(0t/lAt)  =  diag{ai(0t/’/'t) /a2(^t/^^'t)  /  '  •  • }  •  (5.9) 

The  set  :n=l  ,2 ,  •  •  •  ,12}  represents  the  components  of  the 

array  steering  vector  a(0t/’At)  •  Vector  t  (0)  is  given  by 

t(0)  =  (exp[j0J  ,exp[j02]  /  •  ‘  ' /exp[j0i2]  )’'  (5.10) 

and  is  constructed  with  the  element  phases . 

The  rearranged  expression  (5.8)  makes  the  array  steering 
phases  in  a(6t/l/'t)  behave  like  element  phases.  It  also  makes  t(0) 
behave  like  an  array  steering  vector. 

Let  a  transformation  matrix,  dependent  on  the  calibration 
source  direction  (0t;/0t)/  be  defined  as 

T(g,d,,rP^)  A  G'Mg)  .  (5.11) 

This  matrix  has  the  property  that 

r(g,0t/)^t)  b(g, 0,0,0)  =  A'^(0t/0t)  A(0,0)  t(0).  (5.12) 

The  vector  on  the  right-hand  side  has  components  with  unit 
magnitudes .  It  can  be  interpreted  as  an  array  steering  vector  for 
a  new  and  usually  non-physical  'direction'  associated  with  the 
transformation.  For  ( 0 , 0)  =  ( 0^/ 0t)  /  this  vector  becomes  t  (0)  and  is 
independent  of  (0t,0t). 

Applied  to  the  calibration  data,  the  transformation  converts 
each  array  snapshot  to  a  new  one  given  by 
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y  =  r(g,0t/^t)  X 


=  r(g,0t,i/'t)  tee  b(g,<#),  0c/^t)  +  interference  +  noise] 

=  a  A{(9t,^t)  t{0)  +  z 

=  a  T (0)  +  z  (5.13) 

^]^02-g  2  represents  the  contributions  of  the  interference  and.  noise . 
The  first  terra  is  the  transform  of  the  direct  signal  from  the 
transmitter.  This  signal  has  an  amplitude  a  and  a  non-physical 
fixed  'direction'  determined  by  t(0).  The  second  term  is  not 
correlated  with  the  first,  because  the  direct  signal  is  not 
correlated  with  the  interfering  signals  in  the  original  snapshots. 
Moreover,  the  direction  of  each  interfering  signal  in  this  term 
depends  on  as  well  as  the  original  direction  of  the 

interfering  signal . 

Because  the  first  terra  on  the  right -hand- side  of  (5.13)  has  a 
fixed  direction  determined  by  t (0) ,  one  can  used  a  direction 
finding  method  to  estimate  the  element  phases  from  the  transformed 
array  snapshots .  A  technique  that  allows  the  user  to  modify  the 
individual  phases  in  the  array  steering  vectors  must  be  used.  The 
'directional  power  spectrum'  it  produces  is  a  function  of  the  trial 
element  phases . 

5.3.3  Estimation  of  Element  Phases 

The  Capon  method  [9]  is  used  in  direction- finding .  From 
(5.13),  the  data  covariance  matrix  in  this  calculation  is 

Ry  =  E{y  y“} 

=  \a\^  t(0)  +  K,  (5.14) 
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where  E{-}  denotes  the  expected  value,  | a | ^=E{  | q;„| ^} ,  and  R5.=E{zz*^}  . 
The  second  line  strongly  indicates  that  a  directional  spectrum 
calculated  with  Ry  would  have  a  strong  peak  in  the  direction 
corresponding  to  t (0) .  A  more  detailed  analysis  is  given  below. 

Let  0  be  a  trial  value  of  0  in  the  calculation.  The  output 
power  calculated  with  the  Capon  method  is 

1 

(P(0,Ry)  A  -  ,  (5.15) 

F(0,Ry,0) 

where 

F(0„Ry,0b)  =  T»(0J  RyW{0b).  (5.16) 

Because  Ry  is  given  by  (5.13),  application  of  the  matrix  inversion 
lemma  leads  to 


R'/  T  (0)  t”(0)  R;^ 

R-y^  =  r;^ - ,  (5.17) 

1  +  |a|2  F(0,R^,0) 

and 

1  +  |q;|2  F(0,R2,0) 

(P(0,Ry)  =  -  ,  (5.18) 

F(0,R,,0)  +  \a\^  H(0,R,,0) 


with 


H(0,R,,0)  =  F(0,R,,0)  F(0,R,,0)  -  |F(0,R,,0)1" 

=  llqll^*  ||q||^  -  |q“  qP  .  (5.19) 

Here,  ||  •  ||  represents  the  Euclidean  norm.  The  second  line  is 
obtained  by  noting  that  R"^  is  the  inverse  of  an  NxN  non- singular 
correlation  matrix  when  noise  is  present  and  that  it  has  the 
decomposition 
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r  ;1  =  V“  V, 


(5.20) 


3_3  0  nori” singu.l3.ir  NxN  tnstnx  3lso.  This  givsS/  foir  sxsinpls, 

=  t^{4>)  V”  V  t(0) 

=  q«  q  ,  (5.21) 


where 


q  =  V  t(0) 


(5.22) 


and 


q  =  V  t(0)  .  (5.23) 

From  (5.11)  and  the  Schwartz's  inequality 

H(0,R,,$)  £  0.  (5.24) 

The  equality  sign  holds  if  and  only  if  $=<^- 

The  above  results  indicate  that  one  can  obtain  an  estimate  of 
<j>  by  maximizing  (P(0,Ry)  with  respect  to  0,  provided  that  the  rate 
at  which  changes  in  the  neighbourhood  of  <l>  is  small 

compared  with  that  of  |apH(0,R3,0)  .  Increasing  lal  yields  a  more 
accurate  estimate  of  <f>.  Making  F(0,R2,</>)  independent  of  <f>  leads  to 
0=0  at  the  maximum. 

5 . 4  Procedure 


The  procedure  to  calculate  the  element  gains  and  phases  has 
five  steps.  Step  5  involves  the  maximization  of  an  objective 
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function,  and  the  implementation  of  this  step  is  discussed 
separately. 

5.4.1  Steps  in  Procedure 

Let  {x„:  m=l,2,-*-,M}  be  a  set  of  M  array  snapshots.  Let  the 
corresponding  transmitter  directions  be  {  (0n,,^ra)  :  m=l,2,  •  •  •  ,  M}  .  The 
procedure  to  estimate  the  element  gains  and  phases  is  given  below. 

Step  1 . 

Calculate  the  estimates  of  the  element  gains  with  (5.7). 
Denote  the  results  by  {§„:  n=l , 2 , • • • , 12 }  . 

Step  2 . 

Transform  the  snapshots  to  a  new  set  {y„:  m=l,2,--*,M}  with 

Yra  ~  ^  >  m=l ,2,‘‘",M,  (5.25) 

where  T{g,6^,\pJ  is  given  by  (5.11),  and  g  denotes  the  set  of 
gain  estimates. 

Step  3 . 

Calculate  an  estimate  of  the  data  covariance  matrix  as 

1  M 

K  =  -  I  Yn.  yS  (5.26) 

M  m=l 


Step  4 . 

Define  a  function  (P(0,^)  with  (5.15)  and  (5.16),  using 
instead  of  Ry. 
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step  5 . 


Identify  the  estimate  $  as  the  0  that  maximizes  . 

5.4.2  Implementation  of  Step  5 

Step  5  in  the  procedure  requires  the  location  of  the  global 
maximum  of  (P(<f>, 'Ry)  .  In  a  computer  simulation  to  study  the  problems 
that  may  arise,  the  following  results  are  observed: 

1.  The  directional  power  (P(0,Ry)  ,  as  a  function  of  0,  is  usually 
not  unimodal.  However,  the  global  maximum  is  usually  the  one 
closest  to  the  true  (f). 

2.  The  finite  precision  in  computer  arithmetic  adds  a  small 
random  error  to  the  true  value  of  (P(0,Ry)  .  This  error  is  very 
difficult  to  handle  in  the  search  for  the  global  maximum, 
unless  the  constraint  (5.2)  is  imposed  on  every  0  in  the 
calculation  of  the  function.  The  random  error  also  generates 
many  false  maxima  in  the  neighbourhood  of  each  true  maximum. 
The  density  and  peak  values  of  these  false  maxima  generally 
inojfsase  as  the  Euclidean  distance  from  the  true  maximum 
decreases  or  as  the  value  of  the  true  maximum  increases. 

Based  on  the  above  observations,  the  search  for  the  global 
maximum  is  divided  into  two  steps  as  follows: 


Step  5A. 


Initialize  the  components  of  ^  as  samples  of  a  random  variate 
uniformly  distributed  in  the  range  [-90®,90°] ,  assuming  that 
the  true  element  phases  do  not  deviate  from  zero  by  more  than 
90° .  Use  hill -climbing  with  random  searches  at  each  local 
maximum.  Modify  the  components  of  $  individually.  Impose  (5.2) 
on  each  modified  <j)  before  calculating  the  corresponding 
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(P(0,Ry)  . 


Step  5B. 

Repeat  Step  5A  K  times,  where  K  is  a  large  number.  Denote  the 
maximum  positions  and  power  values  by  {  (^k/^Pk)  =  k=l,2,  •  •  •  ,K}, 
respectively.  Identify  the  largest  as  the  value  at  the 
global  maximum  and  denote  it  by  <?.  Identify  the  corresponding 
0k  as  the  estimate 
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6.0  GAIN  AND  PHASE  ESTIMATES 


A  subset  of  approximately  three  thousand  array  snapshots  is 
extracted  from  each  preprocessed  file  produced  in  Chapter  4 .  The 
element  gains  and  phases  obtained  from  these  subsets  are  then 
studied. 

The  procedure  to  construct  a  subset  initially  isolates  a  file 
section  that  spans  360°  in  bearing.  Next,  it  generates  a  set  of  M 
dummy  bearings  as 

=  01  +  360(m-l)/M  .  m=l, 2, • • • ,M-1  ,  (6.1) 

where  6^  is  the  bearing  of  the  first  snapshot  in  this  isolated 
section.  Finally,  it  identifies  the  snapshot  with  bearing  closest 
to  6^  as  x„  and  denotes  the  direction  of  this  snapshot  as  {0m/ ’/'mi¬ 
lt  has  been  stated  in  Chapter  4  that  three  corrections  have 
been  ignored  in  data  preprocessing;  ignoring  the  first  two 
corrections  produce  perturbations  in  the  gain  and  phase  response 
patterns  of  the  array  elements;  these  perturbations  are  usually 
larger  if  the  transmitter  elevation  is  higher;  ignoring  the  last 
correction  produces  amplitude  and  phase  errors  at  the  outputs  of 
the  array  elements;  and  these  errors  are  larger  if  the  transmitter 
range  is  shorter.  Consequently,  the  element  gain  and  phase 
estimates  obtained  from  any  file  are  actually  the  average  of  the 
perturbed  gains  and  phases  in  the  file.  Besides,  the  errors  in 
estimating  these  average  values  are  larger  if  the  transmitter 
elevation  is  higher.  In  the  discussion  of  results,  we  therefore 
give  more  credibility  to  the  estimates  calculated  with  Files  1,  2 
and  3,  where  the  nominal  transmitter  ranges  in  Table  2  are  10  km, 
7  km  and  7  km,  respectively,  and  the  corresponding  transmitter 
elevations  are  4.5°,  11°  and  11°,  respectively.  We  are  also  non- 
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committal  in  the  remarks  on  the  estimates  calculated  with  File  6 
where  the  nominal  range  and  elevations  are  2  km  and  49°, 
respect i ve ly . 

In  Tables  3  and  4  are  the  element  gain  and  phase  estimates, 
respectively,  obtained  with  Files  1  to  6 .  Table  5  contains  the 
transmitter  range,  transmitter  elevation,  the  values  of  M  used  to 
calculate  the  estimates,  the  average  of  the  gain  values  in  the  set 
{4/4/ •  •  • /a2}  /  the  difference  between  the  largest  and  the 

smallest  values  in  this  set.  Some  of  the  more  important 
observations  and  remarks  are  given  below. 

1.  The  standard  deviations  of  the  errors  in  the  gain  estimates 
from  File  2  and  File  3  are  approximately  0.004.  Those  of  the 
phase  estimates  are  approximately  0.2°. 

Remark:  These  standard  deviations  are  derived  from  the 

differences  in  the  estimates  obtained  from  the  two  files.  The 
differences  come  from  the  measurement  errors  in  the  estimates, 
because  the  transmitter  ranges  and  elevations  are  the  same  in 
these  files. 

2.  The  elevation  gain  patterns  of  the  elements  are  different. 

Remark:  One  evidence  of  this  property  comes  from  a  direct 

comparison  of  the  gain  estimates  in  Table  3 .  Another  evidence 
comes  from  Table  5,  where  ave{g,}  is  the  average  of  the  gain 
values,  diff{gi}  is  the  difference  between  the  largest  and  the 
smallest  gains,  and  diff/ave  is  the  ratio  diff  {g,} /ave(gi} .  The 
ratio  diff/ave  increases  steadily  from  0.190  for  File  1  to 
0.255  for  File  6. 

Remark:  This  difference  may  indicate  a  real  dependence  of 

relative  element  gains  on  elevation. 
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3 .  The  element  phases  are  essentially  the  same  in  the  first  five 
files . 

Remark:  This  observation  is  based  on  a  study  of  the  phase 

estimates  in  Table  4 .  The  largest  and  the  smallest  estimates 
in  each  row  differ  by  3.5°  or  less. 

4.  The  phase  estimates  for  elements  n=10  to  12  are  significantly 
different  from  those  of  the  other  elements. 
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Gain  Estimate 

n 

File  1 

2 

3 

4 

5 

6 

(i/'=4 .5°) 

(11°) 

(11°) 

(18°) 

(31°) 

(49°) 

1 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

2 

0.9770 

0.9948 

0 . 9948 

0.9757 

0.9743 

0.9852 

3 

1.0109 

1.0184 

1.0195 

0.9949 

1.0020 

1.0295 

4 

1.0049 

1.0122 

1.0184 

0.9969 

1.0071 

1.0348 

5 

1.1005 

1.1330 

1.1428 

1.1069 

1.1108 

1.1477 

6 

1.0189 

1.0501 

1.0596 

1.0201 

1.0225 

1.0665 

7 

0.9581 

0 . 9938 

1.0010 

0 . 8970 

0 . 9137 

0 . 9198 

8 

1.0119 

1.0481 

1.0513 

1.0262 

1.0410 

1.0496 

9 

0.9113 

1.0706 

1.0729 

1.0988 

1.1519 

1.1867 

10 

0.9980 

0.9836 

0.9876 

0.9828 

0.9866 

1.0169 

11 

0.9790 

0.9979 

0.9979 

1.0292 

1.0924 

1.1350 

12 

0.9432 

0.9365 

0.9434 

0.9223 

0.9373 

0.9926 

Table  3.  Gain  estimates  derived  from  the  first  six  5.1 
MHz  files.  The  target  elevations  are  in  brackets. 
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Table  4.  Phase  estimates  derived  from  the  first  six  5.1 
MHz  files. 
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File 

1 

2 

3 

4 

5 

6 

Range (km) 

10 

7 

7 

5 

3 

2 

'/'(deg) 

4.5 

11 

11 

18 

31 

49 

M 

3600 

3600 

3600 

3600 

3600 

2435 

ave{4} 

0.993 

1.020 

1.024 

1.004 

1.020 

1.047 

diff  {4} 

0.189 

0.196 

0.199 

0.210 

0.238 

0.267 

dif f /ave 

0.190 

0.192 

0.194 

0.209 

0.233 

0.255 

Table  5.  Range,  true  target  elevation  number  of  snapshots  used 
in  estimation  M,  the  average  of  the  gain  estimates  in  the  set 
{4./^/ •  •  •  / a.2}  /  the  difference  between  the  largest  and  the  smallest 
values  in  this  set,  and  the  difference  divided  by  the  average. 
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7.0  EVALUATION  OF  CALIBRATION  METHOD 


The  gain  and  phase  estimates  obtained  from  File  2  are  used  to 
calibrate  the  preprocessed  snapshots  in  Files  1  to  6 .  The  effect  of 
calibration  on  the  performance  of  the  antenna  array  is  then 
studied.  This  study  includes  the  accuracy  of  direction  estimates, 
the  signal-to-interference-plus-noise  ratio  (SINR)  in  the  array 
snapshots,  and  the  resolution  of  signal  directions. 

The  file  subsets  constructed  in  Section  6  are  used  in  this 
study.  They  have  the  convenient  properties  that  each  subset  has 
nsither  too  many  nor  too  few  snapshots;  the  transmitter  circles  the 
array  once;  and  the  snapshots  are  approximately  equally  spaced  in 
bearing. 

In  the  convention  used,  the  calibrated  snapshot  derived  from 
an  uncalibrated  snapshot  x  is  denoted  by  x°.  Its  components  are 
denoted  by  {x?, x^,  •  •  • , and  are  calculated  as 

x=  =  g'^  exp[-j$n]  ,  n=l,2,---,12  .  (7.1) 

Here,  {gi,g2/ *  *  * /gi2}  and  ‘  ’  >^12}  are  the  gain  and  phase 
estimates,  respectively,  obtained  from  File  2.  Their  values  are 
given  by  Column  3  in  Tables  3  and  4,  respectively. 

Although  File  6  is  included  in  the  evaluation,  we  avoid  making 
many  comments  on  the  results  obtained  with  this  file,  because  the 
transmitter  range  is  too  short  and  the  transmitter  elevation,  too 
high. 


7.1  Effect  of  Calibration  on  Direction  Estimates 

This  study  assumes  that  there  is  only  one  signal  in  each  array 
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snapshot  and  that  this  signal  is  the  direct  signal  from  the 
transmitter.  It  uses  the  maximum  likelihood  method  to  estimate  the 
bearing  and  elevation  of  this  signal  in  the  snapshots . 

In  the  presentation  of  results,  the  error  in  a  bearing 
estimate  is  defined  as 

e{d)  =  d  -  e  ,  (7.2) 

where  d  is  the  estimate  and  6  is  the  GPS  bearing  calculated  with 
(3.2)  and  (3.3)  . 

The  errors  in  bearing  estimates  are  plotted  in  point  mode  in 
Figures  13  to  18 .  Their  mean  and  standard  deviations  are  given  in 
Table  6.  Figures  19  to  24  contain  the  elevation  estimates.  Some 
notable  observations  are  given  below. 

1.  The  estimates  seem  to  fall  on  a  smooth  curve,  although  they 
are  plotted  in  point  mode. 

Remark:  It  appears  that  the  SINR  in  the  preprocessed 

snapshots  is  high  and  that  the  interfering  environment  does 
not  change  rapidly  with  bearing. 

2.  Before  calibration,  the  errors  in  bearing  estimates  depend 
strongly  on  the  bearing.  In  particular,  the  errors  are  above 
the  average  value  near  the  60°  bearing  and  are  below  it  near 
240°.  After  calibration,  only  traces  of  this  dependence 
remain . 

3 .  The  mean  value  of  the  bearing  errors  from  each  file  is  not 
zero.  It  is  also  file  dependent. 

Remark:  It  appears  the  non- zero  mean  values  originate  from 

the  assignment  of  A'  in  (3.3)  .  One  can  make  them  equal  to  zero 
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by  using  the  new  values,  denoted  by  A' (new),  in  Table  6. 

Remark:  The  dependence  of  A'  (new)  on  file  number  has  no 

statistical  significance.  After  calibration,  for  example,  the 
difference  between  the  largest  and  the  smallest  values  of 
A' (new)  is  0.7  second  (=  difference  of  values  from  Files  3  and 
6)  .  This  difference  corresponds  to  bearing  errors  of  0.06°  for 
File  3  and  0.15°  for  File  6.  These  errors  are  small  compared 
with  the  corresponding  standard  deviations  in  Column  5. 

4.  The  gross  features  in  Figures  14  and  15  are  very  similar,  if 
the  difference  in  mean  values  is  ignored. 

Remark:  This  property  is  consistent  with  the  interpretation 

that  the  multipath  interference  in  Files  2  and  3  is 
essentially  the  same  and  that  the  variations  in  bearing  errors 
come  mainly  from  the  dependence  of  this  interference  on 
direction. 

5 .  Before  calibration,  the  standard  deviation  of  the  errors  in 
bearing  estimates  are  between  1.09°  and  1.61°.  After 
calibration,  they  are  0.85°  or  less. 

Remark:  The  standard  deviation  after  calibration  appears  to 

be  highly  correlated  with  the  track  deformations  in  Figures  3 
to  8 .  This  standard  deviation  is  the  smallest  in  the  estimates 
from  File  1,  where  the  track  is  a  near  perfect  circle.  It  is 
the  largest  from  File  6,  where  the  track  is  very  much 
distorted. 

6.  The  magnitudes  of  the  elevation  errors  are  generally  smaller 
after  calibration. 

7.  There  are  strings  of  near- zero  elevation  estimates  in  Figures 
19  to  21. 
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Remark:  Inadequate  calibration  may  be  a  reason  for  this 

characteristic.  One  can  verify  by  simulation  that  an  increase 
in  calibration  error  is  usually  accompanied  by  a  reduction  in 
the  magnitudes  of  elevation  estimates,  since  elevation  errors 
have  a  proportional  effect  on  the  elevation  cosine  values 
rather  than  on  the  elevations  themselves. 

8.  The  errors  in  the  elevation  estimates  depend  on  the  true 
elevation.  Usually,  their  magnitudes  are  larger  if  the  true 
elevation  is  lower. 

Remark:  One  reason  for  this  property  is  the  dependence  of  the 
array  elevation  beamwidth  on  elevation.  This  beamwidth  is 
larger  if  the  elevation  is  smaller. 

9.  There  are  more  negative  elevation  errors  than  positive  ones. 

Remark:  One  reason  is  given  in  the  remark  in  Observation  8. 

10.  The  gross  features  in  Figures  20  and  21  are  very  similar. 

Remark:  This  property  is  consistent  with  Observation  4. 
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File 

Mean{ 

(de 

Uncal . 

e{^)} 

ig) 

Cal. 

Std{  e 
(de 

Uncal . 

■(^)} 

jg) 

Cal. 

A'  (n 

(sec< 

Uncal . 

Lew) 

Dnd) 

Cal. 

1 

-0.29 

-0.30 

1.09 

0.32 

7.4 

7.4 

2 

O 

1 

-0.43 

1.19 

0.40 

7.4 

3 

0.61 

0.62 

1.16 

0.39 

7.0 

6.9 

4 

-0.66 

-0.65 

1.23 

0.38 

B 

7.3 

5 

-1.14 

-1.13 

1.29 

0.39 

B 

7.2 

6 

-1.93 

-1.67 

1.61 

0.85 

B 

7.6 

Table  6 .  The  means  and  standard  deviations  of  the  errors 
in  bearing  estimates,  and  the  time  adjustments  needed  to 
make  the  mean  values  equal  to  zero. 
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Figure  13.  Errors  in  bearing  estimates  obtained  with  the  data  in 
File  1 .  The  upper  figure  is  calculated  with  the  uncalibrated  data 
and  the  lower  figure  is  calculated  with  the  calibrated  data. 


Figure  17 .  Errors  in  bearing  estimates  obtained  with  the  data  in 
File  5. 
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Figure  18.  Errors  in  bearing  estimates  obtained  with  the  data  in 
File  6. 
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Figure  19 .  Elevation  estimates  obtained  with  the  data  in  File  1 
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Elevation  (Deg)  Elevation  (Deg) 


Bearing  (Deg) 

Figure  20.  Elevation  estimates  obtained  with  the  data  in  File  2 
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Figure  21.  Elevation  estimates  obtained  with  the  data  in  File  3 
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Figure 


22.  Elevation  estimates  obtained  with  the  data  in  File  4 
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Elevation  (Deg)  Elevation  (Deg) 


Figure  23 .  Elevation  estimates  obtained  with  the  data  in  File  5 . 
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Elevation  (Deg)  Elevation  (Deg) 


Figure  24.  Elevation  estimates  obtained  with  the  data  in  File  6 
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7.2  Effect  of  Calibration  on  SINR 

The  calibration  changes  the  SINR  of  transmitter  signal  in  the 
array  snapshots.  This  change  reflects  the  effectiveness  of  the 
calibration  method,  because  an  increase  in  the  SINR  corresponds  to 
a  higher  detection  probability  of  the  transmitter  signal  and  a  more 
accurate  estimate  of  the  signal  direction.  Here,  the  change  is 
denoted  by  SINRc  and  is  defined  as  the  ratio 

SINRa 

SINRc  =  -  ,  (7.3) 

SINRb 

where  SINR^  is  the  SINR  after  calibration  and  SINRg  is  the  SINR 
before  calibration. 

We  calculate  SINRg,  the  SINR  before  calibration,  with  an 
uncalibrated  array  snapshot  x  as 

SINRb  =  |x///xj" 

=  -  ,  (7.4) 

|xp|a(gt/^t)  P  -  |x«  P 

where  is  the  estimate  of  the  transmitter  direction  obtained 

from  X,  a(@t/^'t)  is  the  array  steering  vector  for  this  direction, 
x.f  is  the  component  of  x  parallel  to  ,  and  Xj^  is  the 

component  of  x  perpendicular  to  .  The  corresponding 

expression  for  SINR^  is  obtained  by  using  x°  and  instead  of 

X  and  t  respectively.  The  justification  for  using  the 

direction  estimates  is  given  in  the  Appendix.  Note  that  the  noise 
power  in  SINRg  and  SINR^  refers  to  the  noise  power  in  the  FFT  bin 
where  the  largest  signal  is  present.  This  bin  has  a  bandwidth  of 
10.4  Hz. 

Figures  25  to  3  0  contain  the  values  of  SINR*,  SINRg  and  SINRc 
calculated  with  the  six  files.  For  the  sake  of  convenience  in  the 
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visual  study  of  results,  a  horizontal  0-dB  dashed  line  has  been 
added  to  the  figures  for  SINRc-  Some  notable  observations  are  given 
below: 

1.  The  average  value  of  SINRg,  the  SINR  before  calibration,  is 

between  16  and  17  dB  in  Files  1  to  5 .  It  is  approximately  15 

dB  in  File  6 . 

Remark:  The  smaller  average  value  in  File  6  is  expected, 

because  the  transmitter  elevation  is  the  highest  in  this  file 
and  the  element  patterns  are  approximately  dipoles  with  nulls 
in  the  vertical  direction. 

2 .  The  average  value  of  SINR^.,  the  change  in  SINR,  is  between  3 

and  6  dB.  The  largest  increases  come  from  Files  2  and  3. 

Remark:  This  property  agrees  with  the  gain  and  phase 

estimates  in  Tables  3  and  4.  The  estimates  obtained  from  File 
2  are  used  to  calibrate  the  array  snapshots.  They  are 
essentially  the  same  as  the  estimates  from  File  3.  Therefore, 
one  expects  the  largest  increase  in  SINR  to  come  from  Files  2 
and  3 . 

3.  The  SINRs  in  Figures  26  and  27,  constructed  with  Files  2  and 
3,  are  very  similar. 

Remark:  This  property  is  also  consistent  with  Observations  4 

and  10  in  Section  7.1. 

4.  In  some  figures,  there  are  directions  where  the  values  of 
SINRc  are  slightly  less  than  0  dB,  thus  indicating  that  SINR^ 
is  slightly  less  than  SINRg. 

Remark:  This  result  is  not  unexpected.  The  values  of  SINR^  in 
these  directions  are  not  unusually  low  compared  with  those  in 
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other  directions.  The  negative  values  only  indicate  that 
original  calibration  coefficients  are  better  for  these 
directions . 

Remark:  The  exceptions  show  that  the  calibration  coefficients 
depend  on  signal  directions. 
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SINRc(dB)  SINRA(dB)  SlNRB(dB) 


Figure  25.  SINRs  calculated  before  and  after  calibration  and 
changes  due  to  calibration,  obtained  with  the  data  in  File  1. 
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Figure  26.  SINRs  calculated  before  and  after  calibration  and 
changes  due  to  calibration,  obtained  with  the  data  in  File  2 . 
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Figure  27 .  SINRs  calculated  before  and  after  calibration  and 
changes  due  to  calibration,  obtained  with  the  data  in  File  3 . 
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Figure  28.  SINRs  calculated  before  and  after  calibration  and 
changes  due  to  calibration,  obtained  with  the  data  in  File  4. 
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Figure  29.  SINRs  calculated  before  and  after  calibration  and 
changes  due  to  calibration,  obtained  with  the  data  in  File  5. 
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7.3  Resolution  of  Signal  Directions 


This  study  compares  the  resolution  performance  of  the  antenna 
array  before  and  after  calibration.  It  involves  the  construction  of 
array  snapshots  with  two  signal  sources  present.  The  maximum 
likelihood  method  is  used  to  estimate  the  two  signal  directions  in 
these  snapshots. 


There  are  five  steps  in  the  procedure.  Steps  1  extracts  a 
section  of  2K  preprocessed  array  snapshots,  , 

from  a  file  subset  constructed  in  Section  6.  This  section  is  60° 
wide  in  bearing  and  is  such  that  the  transmitter  bearing  is  between 
(e"-30°)  and  (0"+3O°),  where  0"  is  the  centre  of  the  region  where 
bearing  resolution  is  studied.  Step  2  separates  this  section  into 
two  groups:  {x'i,x'^,x'^,  •  •  •  ,x'^k.i}  and  {x'^K/*:'k-2/X'k-4/  '  *  ‘  >^"2)  • 
first  group,  the  transmitter  bearing  increases  from  (0"-3O°)  to 
(0"+3O°) .  In  the  second,  the  order  of  the  snapshots  is  reversed  to 
make  the  bearing  decrease  from  (0"+3O°)  to  {0"-3O°).  These  groups 
are  then  combined  to  form  a  new  set  {yi/y2/y3f  '  '  ‘  'yx}  /  which  there 
are  two  signal  sources  crossing  each  other  at  the  bearing  0".  In 
Step  4,  (7.6)  is  obtained  from  (7.1)  by  using  y^  y^  instead  of 

x„  and  x=,  respectively.  The  details  are  given  below: 


Step  1 . 

Extracts  a  section  of  2K  snapshots  from  a  specified  file 
subset.  Identify  its  snapshots  as  {x'l,x'^,x'^,  •  •  •  ,x'^k}  . 


Step  2 . 


Construct  a  set  of  K  new  snapshots  {yi,y2/y3'  ‘  '  '  'Yk) 


yk  =  xVi  +  x" 


2K+2-2k  > 


k=l,2,3, . . . ,K  . 


(7.5) 


Step  3 . 


63 


Use  the  maximum  likelihood  method  to  estimate  two  signal 
directions  from  each  snapshot  in  {yi,y2iY3>  ‘  ‘  ‘  iYk}  •  Denote  the 
estimates  obtained  from  by  {^k2/’^k2}- 

Step  4, 

Construct  a  set  of  K  calibrated  snapshots  {yi,yl,y|,  •  *  • /Yk}  as 
Ykn  =  exp[-j$J  Yta  /  n=l,2,--*,N  .  (7.6) 


Step  5. 

Use  the  maximum  likelihood  method  to  estimate  two  signal 
directions  from  each  snapshot  in  {yl,yl,yl/ *  '  * /Yk}  •  Denote  the 
estimates  from  yl  by  {0ki/^w}  and  {0k2/^k2}* 

This  study  uses  0"=3O°,  K=300  for  Files  1  to  5,  and  K=170  for 
File  6.  The  value  of  6"  corresponds  to  a  bearing  where  the  azimuth 
aperture  of  the  array  is  the  smallest.  Consequently,  the  bearing 
beamwidth  is  the  largest  and  the  changes  in  bearing  resolution  are 
easiest  to  observe.  The  values  of  K  are  determined  by  the  file 
sizes . 

Figures  31  to  36  contain  the  bearing  estimates  @ki^  ^k2»  ^ki/ 
and  §k2  obtained.  These  estimates  are  plotted  in  point  mode.  Each 
figure  has  two  solid  lines.  They  are  GPS  tracks  calculated  with 
(3.2)  and  A '=7. 4  seconds.  One  GPS  track  is  the  source  track  in  {x'i, 

•  ,xVi}  •  The  other  is  the  source  track  in  {x'^k, x V2/ * 2k-4/ 
■-•,x'^}.  There  are  several  notable  observations: 

1.  In  each  figure,  the  bearing  estimates  cluster  around  two 
imaginary  lines  that  intersect  at  bearing  0". 

2.  The  imaginary  lines  deviate  from  the  solid  lines.  The 
calibration  reduces  these  deviations  in  the  results  from  Files 
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1  to  4. 


3 .  The  scattering  of  the  points  about  the  imaginary  lines  is 
generally  larger  if  the  angular  separation  of  the  true  signal 
bearings  is  smaller. 

Remark:  This  property  is  expected,  because  it  is  more 

difficult  to  resolve  two  signal  directions  if  these  directions 
closer  togther.  More  difficulty  in  resolution  also  leads 
to  larger  errors  in  the  estimates. 

4.  The  calibration  reduces  the  scattering  of  the  bearing 
estimates  obtained  with  Files  1  to  5 . 

We  also  guantify  resolution  performance  with  the  technique 
described  below.  It  should  be  noted  that  SINR3  is  at  least  15  dB  in 
each  snapshot  in  the  set  •  •  •  ,xVi}  (Observation  1,  page 

55)  ;  the  SINR;^  in  the  each  calibrated  snapshots  is  higher 
(Observation  2,  page  55);  the  array  aperture  for  0"=3O°  and  zero 
elevation  is  approximately  113  meters;  the  signal  wavelength  at  5.1 
MHz  is  58.8  meters;  and  the  azimuth  beamwidth  of  the  array  is 
approximately  30° .  Therefore,  one  can  say  that  the  signals  in  yj^, 
for  example,  are  not  resolved  if  or  @^2  more  than  15°  from 
and  6^2-  Here,  and  6^^  are  the  GPS  directions  in  x'^^.i  and  x'^K*2-2k/ 
respectively,  and  are  calculated  with  (3.2)  and  A'=7.4  second.  From 
(7.5),  x'^k-i  and  x'^K^.2.2k  are  the  snapshots  used  to  construct  yj,. 

In  the  procedure  below,  and  are  the  number  of  snapshots 
with  unresolved  signal  directions.  The  subscripts  'B'  and  'A' 
indicate  before  and  after  calibration,  respectively. 

Step  1 . 

Let  Jb=0 
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step  2 . 


Do  for  k=l  to  K 

Increase  Jb  by  unity  if 

{|^ki-^ki|  >  15°  and  l^ki-^kzl  >  15°} 
or 

{|^k2-^kil  >  15°  and  I§k2-0k2|  >  15°}. 

End  do 

Step  3 . 

Let  Ja=0 
Step  4 . 

Do  for  k=l  to  K 

Increase  J*  by  unity  if 

{|^ki-0kil  >  15°  and  >  15°} 

or 

{|0k2-ekil  >  15°  and  I@°k2-0k2l  >  15°}. 


End  do 

The  values  of  Jb  and  thus  obtained  are  in  Rows  4  and  5  of 
Table  7,  respectively.  Rows  6  contains  the  percentage  reduction 
(PR)  in  the  number  of  unresolved  snapshots,  i.e.. 
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100  (Jb  -Ja) 


(7.7) 


PR 


Jb 


%. 


■phe  following'  aire  some  of  the  obsen'vations  in  addition  to 
those  derived  from  Figures  31  to  36. 

5.  The  calibration  improves  the  resolution  performance  of  the 
antenna  array  in  Files  1  to  5.  There  is  a  deterioration  in 
resolution  performance  in  File  6 . 

6.  The  values  of  PR  for  Files  1  to  3  are  between  56  and  60%. 
Those  for  Files  4  and  5  are  significantly  smaller.  In  the  case 
of  File  6,  PR  is  negative,  because  there  are  more  unresolved 
signals  after  calibration. 

Remark:  This  dependence  of  PR  on  file  is  not  unexpected.  It 

was  noted  earlier  in  this  section  that  there  are  perturbations 
in  the  gain  and  phase  responses  of  the  array  elements  and  that 
these  perturbations  are  larger  if  the  transmitter  range  is 
shorter  or  if  the  transmitter  elevation  is  higher.  The  range 
is  the  shortest  and  the  elevation  is  the  highest  in  File  6. 
Therefore,  the  poor  results  from  File  6  are  not  unexpected. 
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Bearing  Estimates  (deg)  Bearing  Estimates  (deg) 


Figure  31.  Bearing  estimates  (points)  calculated  with  the  data  in 
File  1.  The  solid  lines  represent  the  true  signal  tracks  calculated 
with  (3.2)  and  A' =7. 4  seconds. 
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Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 
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Bearing  Estimates  (Deg)  Bearing  Estimates  (Deg) 


Figure  35.  Bearing  estimates  (points)  calculated  with  the  data  in 
File  5 . 
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(Deg) 


File 

1 

2 

3 

4 

5 

6 

\p  (deg) 

5.2 

12.3 

12.3 

20.1 

32.0 

50.6 

r  (km) 

10 

7 

7 

5 

3 

2 

Jb 

25 

25 

28 

25 

27 

18 

Ja 

10 

11 

12 

23 

18 

24 

PR  (%) 

60 

56 

57 

8 

33 

-33 

Table  7 .  Number  of  snapshots  with  unresolved  signal 
directions,  and  the  percentage  reduction  (PR)  in  the 
number  of  unresolved  snapshots  after  calibration. 


8 . 0  CONCLUSIONS 


The  authors  have  carried  out  an  analysis  of  the  5.1  MHz 
calibration  data  measured  with  a  high-frequency  antenna  array  at 
Leitrim.  The  objective  was  to  study  whether  the  array  calibration 
method  developed  by  Hung  in  [8]  is  adequate  for  calibrating  the  HF 
antenna  arrays  in  the  CFSRS  modernization  project. 

The  following  results  were  observed  in  the  analysis  of  the 
data  files  generated  with  transmitter  ranges  not  less  than  7  km  and 
transmitter  elevations  not  higher  than  11° . 

1.  With  few  exceptions,  the  calibration  method  increased  the 
signal-to-interference-plus-noise  ratio  (SINK)  in  the  array 
snapshots  by  a  factor  of  2  to  4  (3  to  6  dB) . 

2.  The  calibration  method  increased  the  accuracy  of  target 
direction  estimates.  More  specifically,  it  reduced  the 
standard  deviation  of  the  direction  estimates  by  a  factor  of 
3.0  or  larger. 

3 .  The  method  reduced  the  number  of  snapshots  with  unresolvable 
equal -strength  targets  in  Files  1  to  5 .  The  reduction  is  56% 
or  more  for  the  first  three  files. 

It  is  not  possible  to  decide  whether  the  method  in  [8]  is 
adequate  for  calibrating  the  HF  array  antennas  in  the  CFSRS 
modernization  project.  This  decision  requires  a  study  of 
calibration  data  generated  with  transmitter  ranges  not  less  than  7 
km  and  transmitter  elevations  higher  than  11°.  Data  files  with 
these  characteristics  are  not  available  in  this  study. 

The  difficulties  in  using  the  higher  elevation  files  prompted 
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the  authors  to  study  an  alternative  calibration  method  [11]  that 
tolerates  the  dependence  of  calibration  coefficients  on  directions. 
Preliminary  results  indicated  that  this  alternative  method  was 
superior  to  the  method  in  [8]  .  Because  of  these  results,  the 
authors  have  decided  to  stop  further  work  on  the  method  in  [8]  and 
direct  the  remaining  efforts  to  the  evaluation  of  the  method  in 
[11]  . 
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APPENDIX 


JUSTIFICATIONS  FOR  USING  DIRECTION  ESTIMATES  IN  (7.2) 

The  true  transmitter  direction  is  not  used  in  (7.2)  to 
calculate  SINR^,  SINRb  and  SINRc.  This  direction  cannot  be 
determined  precisely  from  the  recorded  GPSl  and  GPS2  times.  These 
times  are  not  equal  and  there  are  uncertainties  in  the  value  of  the 
time  adjustment  A'  one  should  use  in  (3.2)  . 

In  a  separate  study,  the  following  approaches  to  determine  the 
transmitter  direction  were  evaluated: 

Method  1 : 

Identify  the  transmitter  direction  before  calibration  as  the 
GPS  direction  calculated  with  (3.2)  and  A '=8. 7.  Identify  the 
direction  after  calibration  as  the  GPS  direction  calculated 
with  (3.2)  and  A' =7.4. 

Method  2 : 

Identify  the  transmitter  direction  before  and  after 
calibration  as  the  GPS  direction  calculated  with  (3.2)  and 
A'=8.7. 

Method  3 : 

Identify  the 
calibration  as 
A'=7.4. 

Method  4 : 

Identify  the  transmitter  direction  before  calibration  as  the 
direction  estimated  with  the  uncalibrated  array  snapshots. 
Identify  the  direction  after  calibration  as  the  direction 
estimated  with  the  calibrated  array  snapshots. 


transmitter  direction  before  and  after 
the  GPS  direction  calculated  with  (3.2)  and 
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The  decision  to  present  only  the  results  obtained  with  Method 

4  are  based  on  the  following  observations: 

1.  When  the  GPS  direction  is  used,  the  values  of  SINR^^  and  SINRg 
depend  on  A'.  Using  A'=7.4  instead  of  A'=8.7  could  change 
these  values  by  1.0  dB  or  more. 

Remark:  One  can  show  that  the  dependence  of  SINR^  and  SINRg  on 
A'  originates  from  the  dependence  of  the  errors  e{6)  and  e(^) 
on  A'  .  A  larger  error  results  in  a  smaller  value  of  the 
corresponding  SINR. 

2.  With  Method  4,  the  values  of  SINR^  and  SINRg  do  not  depend  on 
A'  . 

3.  With  Method  4,  a  positive  value  of  SINRc  corresponds  to  a 
higher  probability  of  detecting  the  transmitter  signal  after 
calibration. 

4.  The  values  of  SINR*,  SINRg  and  SINR^  calculated  with  the 
Methods  1  to  3  are  usually  within  2  dB  of  those  obtained  with 
Method  4 . 
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